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We predict differential decay distributions for <j> — > pyy and <fi ~ * ^77 using chiral perturbation theory. We also consider the 
isospin violating decay cj> — > ujty . Experimental information on these decays can be used to determine couplings in the heavy 
vector meson chiral Lagrangian. 

ERRATUM: It was shown by P. Ko et al., Phys. Lett. B366 (1996) 287, that there is a rj pole contribution 
that dominates over what we calculated. 



The 4> factory being built at Frascati is expected to 
produce 10 11 0's per year [Q, enabling rare decay modes 
of the cj> to be measured for the first time. Some of these 
rare decays are of the form V — > V'X, where V and V 
are vector mesons, and can be studied using chiral per- 
turbation theory because the V — V mass difference is 
small. The chiral Lagrangian for V — ► V'X transitions 
written in terms of heavy meson fields was given in a pre- 
vious paper ||. This form of the chiral Lagrangian for 
matter fields has a well-defined momentum expansion jj) , 
and can be used as the starting point for a systematic ex- 
pansion in powers of momentum. The vector meson chi- 
ral Lagrangian has coupling constants (analogous to <7a) 
which have not been measured. In this paper, we show 
that one can determine the unknown coupling constants 
using the rare decay modes 4> — > W77 and <f> — > P77. 

The .vector meson fields can be written as a 3 x 3 octet 
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The chiral Lagrange density which describes the inter- 
actions of the vector mesons with the low-jjnomentum tt, 
K and 77 mesons has the general structure! 
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At leading order in the derivative and quark mass expan- 
sions, 
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where is the vector meson four-velocity, 

VO" = fl"©" + [V*\ 0"] , (8) 
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M is the quark mass matrix M = diag (m u , m^, m s ), 
and 



*We will use the same notation as Ref. M. The notation 
for the pseudoscalar meson fields is standard^ and will not be 
repeated here. The pion decay constant is / = 132 MeV. 



'Here we neglect vector meson widths. These can be taken 
into account by including antihermitian terms in C. 



1 



(10) 



The terms in £kin appear with minus signs because the 
polarization vector is spacelike. 

In the large- N c limit, the octet and singlet mesons can 
be combined into a single "nonet" matrix 
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which enters the chiral Lagrangian. The kinetic, interac- 
tion and mass terms at leading order in 1/N C are 

C kin ^-iTrNUvV)N», (12) 



The value of h p alone does not determine the coupling 
constants g\ and gi- h p vanishes in the large- N c limit, so 
the value of h p is sensitive to deviations of By and gi/g2 
from their large-iV c values. In addition, the chiral loop 
correction to the (fipn coupling may be significant j^j . The 
smallness of h p implies that g\ / gi is close to its large- N c 
value of 2/i/3- 

Like the decay D* — > D s ir°, at leading order in chiral 
perturbation theory the amplitude for the isospin violat- 
ing decay <fi — ► lutt arises from a tree diagram involving 
rfTr mixing H] . This Feynman diagram gives 
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and 
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Comparing with Eqs. (||)-(@), one finds that in the N c 
oo limit, 
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the \(j>) state becomes "pure" \ss), and the nonet matrix 



is 



N u = 
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In studying <f> decays, it is important to determine the 
mixing angle By and the couplings g\ and gi. The mix- 
ing angle By can be determined to first order in flavor 
SU(3) symmetry breaking from the measured vector me- 
son masses 01, 



tan By = 



<m<j, 



- x m K , 
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where the positive sign has been chosen for the square- 
root. The mixing angle is close to the so-called magic 
mixing angle tanBv = 1/V2 ~ 0.71, where the 4> meson 
is a pure ss state. 

The measured <f> — ► pit branching ratio gives | h p \ = 
0.05, where 



where the 4>ur] coupling is 
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Using Q (m d — m u )/m s = 1/44 in Eq. ( |20| ) implies a 
branching ratio, Br(0 — > W7r) = 4.4 x 10~ 6 (fe^ /0.05) 2 . In 
principle the value of the angle By in Eq. (18), the mea- 
sured (f>p7T coupling and a measurement of the branching 
ratio for <fr — > oj7t would determine 31,2- However cor- 
rections from higher orders in chiral perturbation theory 
may be important. It is desireable to have an alternate 
determination of these couplings. Some other cj> decays 
which might have been used to determine 51,2, such as 
4> — > K*K, are not kinematically allowed. 

The coupling constants 31,2 can be determined from 
the radiative <f> decays <f> — > P77, and <p — > ^77. At 
leading order in chiral perturbation theory, there are 
two contributions; the loop graphs in Fig. [TJ and the 
anomaly graph of Fig. ^, and there are no counterterms 
at this order. The computation of the loop graphs is 
straightfoward. We will evaluate the graphs in the limit 
that the photon momentum is small compared with the 
meson mass. In the large-iV c limit, the non-zero loop 
graphs have i-T-meson loops, and the photon momentum 
is smaller than Mr- for the <f> — > P77 and <f> — * luc^t/ de- 
cays. The decay amplitude from the loop graph iall 
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* Vector meson states are normalized to unity. To cqj 
to the usual covariant normalization, multiply Eq. (fez 
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and Ci and hi are the polarization and momentum of pho- 
ton i. The constant A is 
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dx ds 
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where the large- N c values Eq. (26) have been used. The 
differential decay distribution from the anomaly graph is 
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for <fi — > W77 . In the large- N c limit, the constants have 
the simpler form 
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Neglecting isospin nonconservation arising from the 
up-down quark mass difference the amplitude from the 
anomaly graph is 
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M ViTI is the pion or eta mass, h is the (f>pn coupling h p 
(Eq. |lS| ) or the <fxdT) coupling h v (Eq. [2l]) and a PtV — 
1, l/v3 is the anomaly coefficient. Since the parameters 
in the chiral Lagrangian are close to their large- iV c values, 
we will evaluate the radiative decay rate in this limit. The 
anomaly graph and the loop graph are of the same order 
in the l/N c expansion. 

The anomaly amplitude is antisymmetric in <-» e*, 
whereas the loop amplitude is symmetric. As a result, 
there is no interference term for the spin averaged decay 
rate, so we will consider the loop and anomaly decay 
rates separately. It is convenient to use the dimensionless 
variables 
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where k\ and k 2 are the four-momenta of the two pho- 
tons, and Ei and E 2 are their energies in the rest frame 
of the 4> meson. The allowed region in the Dalitz plot 
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in the xs plane is < x < [(1 + r — s 

0<s<(l-Vr) 2 - 

The differential decay distribution from the loop graph 

is 



Integrating over x from to [(1 + r — J) z 
gives 
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Total decay rates are obtained by integrating dT/ds 
from s = to s = (1 — \/r) 2 . The branching ratios 
from the loop graph for (f> — > /J77 and <f> — ► W77 are 
5.8 x 10- 9 (. 92 /0.75) 4 and 4.2 x 10^ 9 (.g 2 /0.75) 4 , respec- 
tively. The non-relativistic quark model predicts that 
g 2 = 1, and the somewhat smaller value g 2 ~ 0.75 is sug- 
gested by the chiral quark model ||. These branching 
ratios are quite small, but the radiative decays should be 
observable at the 4> factory. 

The anomaly graph for cj> — > W77 involves the unknown 
coupling h^, which vanishes in the large- N c limit. Since 
hrf and h p both vanish in the large- N c limit we will use 
h n « hp to estimate the importance of the anomaly 
graph relative to the loop graph. With this assump- 
tion for h v , the anomaly graph is much smaller than the 
loop graph over the entire Dalitz plot, and can be ne- 
glected (The ratio of decay widths is r anoma i y /ri oop = 
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5 x Kr 4 (V0.05) 2 (0.75/g 2 ) 4 -) Thus a measurement of 
the <j> — > W77 decay ratel allows one to determine the 
coupling constant 52, and g\ — 2g2/V%- It might even 
be possible to compare the chiral peturbation theory pre- 
diction for the differential decay distribution with exper- 
iment if <?2 is reasonably large, and the (f> factory achieves 
its design luminosity. 

The situation is quite different for (f> — ► P77. The 
branching ratio from the anomaly graph for <j) — > ^77 
is the same as the branching ratio of 4.3% for the process 
<\> — > p°tt°, since the 7r° can be on-shell, and it decays al- 
most exclusively into two photons. Thus the total decay 
rate from the anomaly overwhelms that from the loop 
diagram. However the rate from the anomaly graph is 
concentrated in a band in the xs plane centered around 
s = Ml/rrA = 0.0175. The ratio of the loop to the total 
decay distribution is show in Figs. (HJH). There are re- 
gions in the Dalitz plot where the loop graph is the dom- 
inant decay mechanism. Measuring the radiative decay 
rate in this region should allow one to determine ri oop , 
and thus the coupling constant <?2- 

The rate for <f> — > P77 near the pion-pole is mainly from 
the anomaly graph. A measurement of the rate in this 
region of the Dalitz plot probes the q 2 dependence of the 
7T° — > 77 amplitude, which is theoretically interesting. 
The q 2 dependence arises from higher order terms in the 
meson chiral Lagrangian that contain an e-symbol, and 
also from loop graphs involving the Wess-Zumino term. 
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§For s very near M^/m^, = 0.0175 the <f> — > 1^77 rate is 
dominated by the isospin violating process 4> — > u>n — > CJ77 
which is expected to have a branching ratio around 4 x 10 
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FIG. 1. Loop Graphs for <f> —> P77 and cf> — > ^77 that contribute at leading order in chiral perturbation theory. Graphs with 
a photon attached to a vector meson line vanish in v ■ e = gauge. 
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FIG. 2. The anomaly graph for <f> — > P77 and <f> — > CJ77. 
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FIG. 3. The ratio dT\ oov /dx ds/(dTi 00 p/dx ds + dr anoma i y /cix ds) over the Dalitz plot. The pion pole is at s — 0.0175. 




